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1 Introduction 



Pairing correlation between nucleons is a key ingredient to describe the struc- 
ture of neutron stars and finite nuclei. There are two distinct ways of descrip- 
tion of such finite- density nuclear many-body systems; the non-relativistic and 
the relativistic ones. The latter incorporates the mesons explicitly in addition 
to the nucleons in terms of a field theory. Both describe the basic proper- 
ties such as the saturation with a similar quality in different manners. Irre- 
spective of whether non-relativistic or relativistic, however, various theoretical 
approaches can be classified into two types: One is realistic studies adopting 
phenomenological interactions constructed for finite-density systems from the 
beginning (hereafter we call this the P-type, indicating "phenomenological"), 
as often done in the studies of heavy nuclei. And the other is microscopic 
studies based on bare nucleon-nucleon interactions in free space (the B-type, 
indicating "bare"). In relativistic studies, typical examples of these two types 
as for the particle-hole (p-h) channel are the relativistic mean field (RMF) 
model and the Dirac-Brueckner-Hartree-Fock (DBHF) method, respectively. 
As for the particle-particle (p-p) channel, that is, pairing correlation, a bare 
interaction was used as the lowest-order contribution in the gap equation [1] 
in a study of the B-type [2]. This is thought to be a good approximation at 
least for the ^5*0 channel (Refs.[3,4], for example). The first relativistic study 
of the P-typc of pairing correlation in nuclear matter was done by Kucharck 
and Ring [5]. They adopted a one-boson exchange (OBE) interaction with 
the couphng constants of the RMF model, which we call the RMF interac- 
tion hereafter, aiming at a fully selfconsistent Hartree-Bogoliubov calculation, 
which we call the Pl-type, in the sense that both the p-h and the p-p channel 
interactions are derived from a common Lagrangian. But the resulting pairing 
gaps were about three times larger than those accepted in the non-relativistic 
studies (see Fig.6). The reason can be ascribed to the fact that the couphng 
constants of the RMF model were determined by physics involving only low 
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momenta {k < k^, -^{kp)^ — po denoting the saturation density of symmetric 
nuclear matter), and therefore the adopted OBE interaction is not rehable at 
high momenta. After a five-year blank, some attempts to improve this were 
done [6-8]. But their results were insufficient. 

An alternative way is to adopt another interaction in the p-p channel while 
the single-particle states are still given by the RMF model. We call this the 
P2-type. There are some variations of this. The first one, which we call the 
P2a-type, adopts another phenomenological interaction for the p-p channel. 
Actually, the non-relativistic Gogny force [9] was used combined with the 
single-particle states of the RMF model for finite nuclei in Ref. [10] and sub- 
sequent works, and gave excellent results. The second variation, which we call 
the P2b-type, is to adopt a bare interaction that describes the high-momentum 
part realistically, the Bonn potential, again combined with the single-particle 
states of the RMF model [11] (see also Ref. [12]). If one assumes that the 
RMF model simulates roughly the DBHF calculation, this P2b-type can be 
regarded as simulating the B-type calculation [2] mentioned above. The re- 
sults of these P2a- and P2b- types are very similar at densities p < Pq. This 
supports a statement that the Gogny force resembles a realistic free interac- 
tion in the low-density limit [13] (see also Ref. [4]). But a clear difference can 
be seen at p ~ po- This difi^erence can also be seen in fully non-rclativistic cal- 
culations; compare the results in Refs.[14] and [15], for example. The precise 
origin of this difference has not been understood well. A comparison after tak- 
ing the polarization effects which have been known to be important at finite 
densities [16-18] into account may be necessary (see also Ref. [12]). The third 
variation, which we call the P2b'-type, is to parameterize the p-p channel in- 
teraction in terms of a few scattering parameters. This was actually examined 
by being combined with the DBHF calculation [19], which we call the B'-type 
(see also Ref. [20]). This method is free from model-dependent ambiguities and 
meets the viewpoint of modern effective field theories [21,22] but is applicable 
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only to dilute systems. These classifications are summarized in Table 1. 

Since we are interested in a wide density range where the ^5*0 gap exists and 
would like to respect the selfconsistency in the sense that both the p-h and 
the p-p channel interactions are derived from a common Lagrangian, here we 
construct a phenomenological relativistic nucleon-nucleon interaction based 
on the RMF interaction (the Pl-type) by adjusting to the pairing properties 
given by the RMF+Bonn calculation (the P2b-type). In other words, we aim 
at constructing an interaction similar to the Gogny force in the sense that it 
reproduces the pairing properties given by the bare interactions in spite of the 
fact that it was constructed for the finite-density system from the beginning. 

The contents of this paper are as follows: In sect. 2 we present our method 
of constructing a phenomenological p-p channel interaction for the super- 
fiuid gap equation. This procedure is applied both to the Hartree and to 
the Hartree-Fock model. In sect. 3, we discuss the pairing properties obtained 
in the Hartree approximation. Note here that we adopt the no-sea approx- 
imation throughout this paper and therefore the Relativistic Hartree model 
means the so-called Relativistic Mean Field model. In sect. 4, the Relativistic 
Hartree-Fock-Bogoliubov calculations with or without the modulation of the 
high-momentum interaction are presented. Conclusions are given in sect. 5. 

2 Construction of a relativistic particle-particle channel interac- 



tion 



We start from the ordinary a-uj model Lagrangian density, 




(1) 
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The antisymmetrized matrix element of the RMF interaction V derived from 
this Lagrangian is defined by, 

^;(p, k) = (ps', pV|y |ks, ks) - (ps', pV|y|ks, ks), (2) 



under an instantaneous approximation, with tildes denoting time reversal. 
After a spin average and an angle integration are performed to project out the 
yS'-wave component, its concrete form is given by 



9l i 2(E;E; + M*')-(p' + e + ml) ( {p + kr + mr 
2E;E*,\ ^ Apk ""[(p-ky + rnl^ 



2E;Elpk 

where 



£;* = Vk2 + M*2. (3) 



Our policy of constructing a phenomenological interaction proposed above is to 
introduce a density-independent parameter A so as not to change the Hartree 
part with the momentum transfer q = which determines the single-particle 
energies, respecting that the original parameters of the RMF are density- 
independent. In the Hartree-Fock model, the q 7^ part also contributes 
to the mean field. This will be discussed in sect. 4 in detail. Since the high- 
momentum part of the RMF interaction does not have a firm experimental 
basis as mentioned above, we suppose there is room to modify that part. Need- 
less to say, such a modification should be checked by studying independent 
phenomena, for example, medium-energy heavy-ion collisions. Some adjust- 
ments of the density-dependent 5 force to the ones that give realistic pairing 
properties have already been examined in the non-relativistic studies [13,4]. 
Among them, a fit to a bare interaction in the T = channel was done [4]. 
We aim at a similar procedure in the relativistic model for the T — 1 channel. 
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In the preceding letter [23], the upper bounds of the momentum integration 
in the gap equation, 



A(p) = — \- Jv{p, k) , ^^^^ edk , (4) 



where 



Ek = El + g^{uj''), (5) 



and the nucleon effective mass equation. 



M* = M - / vlk^dk , (6) 



where the spin-isospin factor 7 = 4 and 2 indicate symmetric nuclear matter 
and pure neutron matter, respectively, were cut at a finite value A, as usually 
done in condensed-matter physics, since the gap increases monotonically until 
reaching Kucharck and Ring's value when the model space is enlarged as 
shown in Fig. 3 of Ref.[8], while v{p, k) is left unchanged. We call this method 
the sudden cutoff hereafter. This was done first in Ref.[5] by inspection. We 
proposed a quantitative method to determine A, which is described below, 
and obtained 3.60 fm~^ for the hnear a-u parameter set in Ref.[23]. This value 
almost coincides with their value, about 3.65 fm~^, for the NL2 parameter set. 

In the present paper, we examine smooth cutoffs that weaken the high-momentum 
part; a form factor /(q^), q = p — k, is applied to each nucleon-meson vertex 
in f;(p, k) while the upper bounds of the integrals are conceptually infinity. 
They are replaced numerically by a finite number, 20 fm~^ which has been 
proved to be large enough in Ref.[8]. Since there is no decisive reasoning to 
choose a specific form, we examine four types. 



A^ 

monopole: /(q^) = — - — - 
A^ -f q^ 
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dipole: f{q^)-- 



strong(a): /(q^ 



A2-q2 



A2 + q2 ■ 

strong(b): /(q') = , {i^a,u;). (7) 

Note that the sudden cutoff above was apphed to k, not to q. 

The parameter A is determined so as to minimize the difference in the pairing 
properties from the results of the P2b-type RMF+Bonn calculation. Assuming 
the P2b-type roughly simulates the B-type as mentioned in sect.l, conceptu- 
ally we aim at fitting to the pairing properties given by the fully microscopic 
B-type calculation. Here we adopt the Bonn-B potential because this has a 
moderate property among the available (charge-independent) versions A, B, 
and C [24]. The pair wave function, 



E^^(k)^^(E,-E,^)^ + A^(k), (8) 
is related to the gap at the Fermi surface, 

oo 

A(A;f) = j v{kF, k)(j){k)k^dk , (9) 



and its derivative determines the coherence length [25], 

^ [ J,^ \cP\^k^dk ) ' 



(10) 



which measures the spatial size of the Cooper pairs. These expressions indicate 
that A(A;f) and ^ carry independent information, (f) and respectively, in 
strongly-coupled systems, whereas they are intimately related to each other 
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in weakly-coupled ones. Therefore we search for A which minimizes 



X 




A(A:f)rmf - A(A;f)b, 

A(A;F)Bonn 



lonn 



) 



2 



+ 



( 



CrMF — '^Bonn 
^Bonn 



) 



2 



(11) 



where the subscripts "RMF" and "Bonn" denote the RMF interaction includ- 
ing A and the Bonn-B potential, respectively, while the single-particle states 
are determined by the original RMF model in both cases. 

The actual numerical task is to solve the gap equation (4) and the effective 
mass equation for the nucleon (6). They couple to each other through Eq.(5) 
and 



The parameters of the standard a-u model that we adopt are g^. — 91.64, 
gl ^ 136.2, rua = 550 MeV, = 783 MeV, and M = 939 MeV [26]. N 
in is taken to be 11; k-p = 0.2, 0.3,..., 1.2 fm^^. In the following, the 
results for symmetric nuclear matter are presented. Those for pure neutron 
matter are very similar except that A(/cf) is a little larger due to a larger 
effective mass M* as shown in Fig. 1 (b) of Ref.[23]. Minimizations of give 
the optimal values, A — 7.26, 10.66, and 10.98 fm~^ for the first three types of 
form factor, respectively, as shown in Fig. 1(a). We do not choose an optimal 
A for the strong(b) type because of its pathological A-dependence shown in 
Fig. 1(b). Figure 1(b) shows that the A-dependence of these smooth cutoff cases 
is very mild, except for the strong(b) type, in comparison with the sudden 
cutoff case. The very steep A-dependence in the strong(b) case is due to the 
consecutive depression of the attraction and the repulsion at around A ~ m„ 
and ~ m^, respectively. In the strong(a) case, although another smaller A 
around 3 fm~^ can give similar A{kF), v{kF, k) and consequently A{k) exhibit 
an unphysical staggering at around k ~ m^. Therefore we discard this. In 
addition, although the strong(a) case with A = 10.98 fm~^ gives practically 




(12) 
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the same /cp-dependence of A(A;f) and ^ (in Fig.2 shown later), v{p,k) — 
at A = |q| = |p — k| brings about an unphysical staggering in A(A;) and 
0(A:); this leads to an oscillatory structure in r-space with a period ~ tt/A. 
Therefore wc discard this, too. Form factors with similar A arc also suggested 
in a study of medium-energy heavy- ion collisions [27]. This indicates that the 
present results have a physical meaning. 

3 Pairing properties obtained in the Relativistic Hartree-Bogoliubov 
calculation by using the constructed interaction 

Figure 2 presents the results for A(A;f) and ^ as functions of the Fermi mo- 
mentum /cp, obtained by using the cutoff parameters so determined. Both 
the monopole and the dipole cases reproduce the results from the Bonn-B 
potential very well, as the sudden cutoff case studied in Ref.[23], in a wide 
and physically relevant density range, in the sense that pairing in finite nu- 
clei occurs near the nuclear surface where density is lower than the saturation 
point [28-30] and that the calculated range of /cp almost corresponds to that of 
the inner crust of neutron stars [31]. This is our first conclusion. In the present 
method with a density-independent cutoff parameter, some small deviations 
remain: The overall slight peak shift to higher kp in A(/i;p) and the deviation 
in ^ at the highest kp are brought about by the systematic deviation in the 
critical density where the gap closes, between the calculations adopting bare 
interactions and those adopting phenomenological ones as mentioned above. 
The deviation at the lowest kp is due to the feature that the present model is 
based on the mean-field picture for the finite-density system; this is a differ- 
ent point from the three-parameter fitting in Refs.[13,4]. Actually, in such an 
extremely dilute system, the effective-range approximation for free scattering 
holds well [19]. 

Next we look into the momentum dependence at kp — 0.9 fm~^, where A(/cp) 
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becomes almost maximum. Figure 3(a) shows (f){k). It is evident that the 
monopole form factor gives the result identical to the Bonn-potential case 
as in the sudden cutoff case. This demonstrates clearly the effectiveness of 
the interaction constructed here. Since we confirmed that the results of the 
monopole and the dipole cases coincide within the width of the line, hereafter 
only the monopole case will be shown as a representative. This quantity peaks 
as seen from Eq.(8). The width of the peak represents the reciprocal 
of the coherence length. Equation (8) shows that 0(/c) is composed of A(/c) 
and the quasiparticlc energy Eqp[k). Figure 3 (b) graphs the former. The gaps 
of both the sudden cutoff and the monopole form factor cases are almost iden- 
tical up to A; ~ 2/cf, and deviations are seen only at larger momenta where 
Eqp{k) are large and accordingly pairing is not important. This is not a trivial 
result since the bare interaction is more repulsive than the phenomenological 
ones constructed here even at the momentum region where A(A;) are almost 
identical as shown in Fig. 3 (c). The reason why we compare the constructed 
p-p channel interaction with the bare interaction is as follows: Although evi- 
dently the RMF interaction corresponds to a medium-renormalized one, not 
to a bare one, here we aim at constructing an interaction similar to the Gogny 
force in the sense that it reproduces the pairing properties given by the bare 
interactions in spite of the fact that it was constructed for the finite-density 
system from the beginning. Accordingly we compare them to see the differ- 
ence between an effective and a bare interactions which give similar pairing 
properties. This figure indicates that the difference is roughly /c-independent. 

Here we turn to the dependence on r, the distance between the two nucleons 
that form a Cooper pair, in order to look into the physical contents further. 
The gap equation, before the angle integration that results in Eq. (4), can 
be Fourier-transformed to the local form, A(r) = — v(r)0(r) in r-space in the 
non-relativistic limit [32]. One can see from this expression that, assuming 
A(r) is finite, 0(r) is pushed outwards when v{r) has a repulsive core, as the 
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Brueckner wave functions [3,33]. This is related to an observation that the 
gap equation reduces to a Schrodinger equation for the relative motion of the 
two particle that form a Cooper pair in the limit of 0, i.e., at high 

/c[34,35,4]. This, on the other hand, masks practically the differences in the 
repulsive interactions at short range, in other words, widely spread in /c-space. 
The r-space pair wave functions, 

= mMkr)edk, (13) 



where jo{kr) is a spherical Bessel function, at kp — 0.9 fm"^ are shown in 
Fig. 4(a). Appreciable differences are seen only in the core region as men- 
tioned above. The coherence length, that is a typical spatial scale of pairing 
correlation, is about 6 fm at this kp as shown in Fig. 2 (b); this is almost one 
order of magnitude larger than the size of the core region. Therefore, practi- 
cally we can safely use the p-p channel interaction, including the sudden cutoff 
one, constructed here for the gap equation. Figure 4 (b) shows the correspond- 
ing A(r). The gaps are positive at the outside of the core and negative inside 
in all cases. Note here that the gap equation is invariant with respect to the 
overall sign inversion; we defined as A{kp) > 0. Their depths at the inside 
region reflect the heights of the repulsive core. In the sudden cutoff case, A(r) 
behaves somewhat differently from others due to the lack of high-momentum 
components. 

4 Relativistic Hcirtree-Fock-Bogoliubov calculation with a cutoff 

In the preceding sections, we have shown that we can construct phenomeno- 
logically relativistic p-p channel interactions which give realistic pairing prop- 
erties by introducing a density-independent momentum-cutoff parameter to 
the (no-sea) Relativistic Hartree model. To do this, we have fully utilized the 
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property that only the q = part of the interaction contributes to the Hartree 
mean field. To see the further applicabihty of the method presented above, the 
Relativistic Hartree-Fock (RHF) model, in which the q 7^ part of the inter- 
action also contributes to the p-h channel, should be examined. In the follow- 
ing, we investigate the Relativistic Hartree-Fock-Bogoliubov (RHFB) model 
with a momentum-cutoff form factor, by comparison with the corresponding 
Hartree-Bogoliubov calculation. An RHFB calculation with a sudden cutoff 
was previously done by Guimaraes et al.[6]. They obtained very large pair- 
ing gaps in the no-sea approximation (Fig. 4 of Ref. [6]). Before studying the 
effects of the form factor to modulate smoothly the high-momentum interac- 
tion, we examine the sudden cutoff — this is a straightforward extension of 
our previous calculation in Ref. [23] to the RHFB. 

The RHF model was described in detail in Refs. [26,36,37], for example. The 
difference from the Hartree model is the second term in Fig. 5. This introduces 

the space component of the vector selfenergy, Tl" , and all the Lorentz compo- 
nents of E — the scalar E*, the vector (time) E°, and the vector (space) E^ 
— become momentum dependent: 



In contrast, in the Hartree model, is not present, S'' and S° are momentum 
independent, and TP is fully determined by the input density. Their concrete 
forms are 



E(p) = E^(p)-70E°(p)+7-pE^(p). 



(14) 




1 



/ 







dkk 



M*{k) rl 
E*{k) ^4 
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/ 



dkk[-gle^{p, k) + -g'Mp, k)\ , 

t k* rl 

2 J dkk^-^[-gl^a(p,k) + gl^^(p,k) 







1 



(15) 







with 



eiip,k) 



Ai{p, k) + 2pk 
Ai{p,k)-2pk ' 



<^i{p,k) 




A{p, k) 



2 I 1 2 I 2 

p + k + , 



(16) 



as Eqs. (5.76) - (5.79) in Ref.[26]. The retardation effect is neglected here. 
Note that the momentum dependence of M* is stemming from that of E*: 
M* = M + (see Eq.(6)). When pairing is introduced, Jq^ is replaced by 
vl, and therefore T,^{p), ^^{p), E^(p), and A(p) (Eq.(4)) couple through 
Eqs. (12) and (5). When the momentum space is discretized to n meshes, these 
quantities form a set of 4n-dimensional coupled non-linear equations; this is a 
contrast to the n + 1-dimensional ones — A(p) and a momentum- independent 
M* or E* — in the Hartree case. 

In order to reproduce the saturation, the coupling constants are adjusted to 
g^. = 83.11, g^ = 108.05 with the masses being unchanged [26]. We adopt this 
parameter set for the time being. Other sets will be examined later. Note that 
the pairing contribution to the energy density around the saturation point is 
neghgible. First we compare the gap at the Fermi surface calculated without 
a cutoff in RHFB and RHB, in Fig. 6. We call them the full calculations here- 
after. The calculated gaps in the former are larger than those in the latter, by 
about 4 MeV at the maximum, for example. This difference is brought about 
by that in the couphng constants rather than the Fock effect itself [38] . Actu- 
ally we confirmed that a Hartree calculation with the coupling constants of the 
Hartree-Fock gave almost the same gaps as those given by the Hartree-Fock 
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calculation although such a calculation destroys the saturation completely. 
This result reflects the property that the gap is not sensitive to the detail 
of the single-particle states. If the coupling constants of Ref.[6], gl = 96.392, 
gl, = 129.260 are adopted, calculated gaps become even larger; about 15 MeV 
at the maximum, for example. But it is evident that this is still smaller than 
their calculated value in their Fig. 4. The reason for this difference is not clear. 

Now let us introduce a sudden cutoff in the upper bound of the momentum 
integrations in S(p)'s and A(p) as in the previous Hartree case. The result 
is graphed in Fig. 7 by the long-dashed curve. Although a plateau appears 
around A = 2 fm~^ as in the corresponding case in Fig. 1 (b), the gap value of 
the present plateau is about 4.5 MeV, which is larger than the physical value 
given by the bare interaction, about 2.8 MeV at this density. Applying the 
same procedure as in the previous Hartree case, the obtained optimal cutoff is 
1.26 fm^-*^. This indicates that the attractive interaction alone accounts for the 
physical magnitude of the gap since the A of the plateau corresponds to the k 
at which v{k-F,k) turns from attractive to repulsive [23]. Note that evidently 
the RHFB calculation with this sudden cutoff can not be applied to the case 
with /cp > A. 

Next we proceed to the smooth cutoff — the form factor in the nucleon- 
meson vertices. The monopole and the dipole types are examined here. Their 
A-dependence is included in Fig. 7. The optimal cutoffs are A = 3.40 fm~^ and 
5.02 fm~^ for the monopole and the dipole types, respectively. These values 
are smaller than the corresponding ones in the Hartree cases since the original 
gap values given by the full calculation are larger. Interestingly, however, their 
ratio, A(monopole)/A(dipole), almost coincides with that of the Hartree case. 
The p-p channel interactions including the optimal monopole form factor are 
compared in Fig. 8. The high-momentum repulsive part, in particular at lower 
densities, is strongly suppressed in the RHFB case. This testifies the discussion 
about Fig. 7 above that the optimal cutoffs are smaller. In Fig. 9 we present 
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A{kp) and ^ calculated by adopting the optimal cutoffs. They reproduce the 
values given by the Bonn-B potential to an extent similar to the Hartree case 
or a httle better except A(A;f) of the sudden cutoff case. Again the results of 
the monopole and the dipole form factors coincide with each other within the 
width of the line. 

Although the pairing properties of the RHFB with the optimal form factor 
presented in Fig. 9 are essentially the same as those of the corresponding 
Hartree calculation, an essential feature of the Hartree-Fock model is that the 
form factor can affect the saturation property. Therefore we have to check 
this before concluding the applicability of the present procedure. The energy 
density is graphed as a function of the Fermi momentum in Fig. 10. This 
figure can be compared with Fig. 40 in Ref.[26]. Here the pairing correlation 
energy (both cases) and the cutoff contribution (Hartree-Fock case only) are 
additionally included. Although these two curves almost coincide with each 
other at low densities, the cutoff effect destroys the saturation in the RHFB 
case. This is because suppressing the high-momentum repulsion breaks the 
balance between the attraction and the repulsion, and the latter contribute 
more at higher densities (see Fig. 8(a)). To see this more closely, the cutoff 
dependence of the energy density without the pairing contribution at the sat- 
uration density, fcp = 1.42 fm^^, is shown in Fig. 11. In the large-A limit, the 
energy densities approaches to about -15.69 MeV, which is 0.38 % less bound 
than the original value in Ref.[26], -15.75 MeV, because of the instantaneous 
approximation. This figure shows that 

Ap_p < Ap_h , (17) 

if we call the cutoff which reproduces the pairing properties and which does 
the saturation properties Ap_p and Ap_h, respectively. 

Up to now, we used a parameter set of Serot and Walecka. To complete the 
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discussion, other sets with different characteristics are also examined here. 
One is a set used by Bouyssy et al., gl — 69.62, — 153.81, — 440 MeV, 
= 783 MeV, and M = 938.9 MeV, which gives saturation at = 1.30 
fm~^[37]. A distinct feature of this set is that the plateau appears at very small 
pairing gap (see Fig. 7). This indicates that the repulsive part also contributes 
to the gap as in the Hartree case, and accordingly, Ap_p will be larger than 
in the case of Serot and Walecka's set. Actually, the optimal monopole cutoff 
is 5.22 fm~^, which is larger than 3.40 fm~^ in the previous case. But Fig. 11 
shows that still Ap_p < Ap_h- Another is a set used by Jaminon et al.[36]. 
An interesting feature of this set is that the saturation (at kp — 1.36 fm~^) is 
given by introducing a "weak" form factor [39], 



from the beginning. Note that here is the square of the 4- momentum transfer 
and therefore this contains a retardation effect. The parameters are gl — 93.87, 
gl = 127.55, m„ = 550 MeV, = 782.8 MeV, and A = 7.754 fm'^ Since M 
is not shown explicitly, 939 MeV is adopted. The RHFB calculation with this 
form factor under the instantaneous approximation gives pairing gaps which 
is larger than the physical values; A(/cf) ~ 8.5 MeV at the maximum, for 
example. This means that still Ap_p < Ap_h is necessary to reproduce the 
physical pairing gap. 

In all the calculations above, the cutoff in the form factors is common to a 
and ui. Referring to the Bonn potential and Ref.[27], for example, Acr 7^ A^^ is 
another option. Of course large deviations from A^. = A^^ destroy the satura- 
tion at least when the coupling constants are kept unchanged. We examined 
some cases in a hmited range, 4^ = 0.70 (~ ^) - 1.2, but the results were 




(18) 



negative. 
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5 Conclusions 



We have constructed phenomenologically a particle-particle channel interac- 
tion which suits the gap equation for nuclear matter. This was done by intro- 
ducing a density-independent momentum-cutoff parameter to the one-boson 
exchange interaction derived from the Lagrangian of the RMF model and ad- 
justing it to the pairing properties obtained by the Bonn-B potential. The 
model pairing properties were calculated by using the RMF model in the p-h 
channel and the Bonn-B potential in the p-p channel utilizing the properties 
that the RMF model simulates crudely[^ the G-matrices in the DBHF calcu- 
lation and that the pairing properties are not sensitive to the detail of the 
single-particle states. By this procedure we aimed at constructing an interac- 
tion like the Gogny force in the sense that it reproduces the pairing properties 
given by the bare interactions in spite of the fact that it was constructed for 
the finite-density system from the beginning. The actual determination of the 
optimal cutoff parameter was done by the method proposed in our preceding 
letter for the Hartree mean field plus a sudden cutoff [23]. 

In the present paper, first we applied this method to the (no-sea) Relativistic 
Hartree model plus various types of form factor which modulates the high- 
momentum repulsion that spoils the pairing properties. Among the four types 
examined, the monopole and the dipole form factors exhibit desired properties. 
Close analyses in momentum space and coordinate space have clarified that 
the gap equation involves a mechanism to mask the difference in the short- 
range repulsion between the bare and the in-medium effective interactions 
and that, in the typical spatial scale of pairing phenomena determined by the 
coherence length, practically there are no difference in the pairing properties. 



^ The coupling constants must be density-dependent in order to simulate the G- 
matrices quantitatively. 



17 



Second we performed RHFB calculations with and without a momentum cutoff 
— a sudden cutoff and two types of form factor. The reason why we examined 
the RHFB separately is that the form factor which modulates the q 7^ part of 
the interaction also affects the p-h channel. The same procedure was applied; 
the resulting optimal cutoffs are smaller than those for the Hartree model 
because the original gap values given by the full calculation without a cutoff 
are larger. Interestingly, the ratio A(monopole)/A(dipole) almost coincides 
with that for the Hartree model. Having confirmed that these optimal form 
factors reproduce A(A:f) and ^ with a quality similar to that in the Hartree 
calculation, we looked at the saturation curve of the energy density. Although 
the results of the RHB and the RHFB each with the optimal form factor 
coincide with each other up to /cp ~ 1 fm"^, the result for the latter starts 
to deviate at larger fcp. This indicates that Ap_p < Ap_h is necessary in order 
also to reproduce the saturation simultaneously using the RMF interaction 
both in the p-h and the p-p channels. This holds also for the other parameter 
sets with characteristic features mentioned in the previous section. 

These results support the observation that the high-momentum components 
of the original RMF interaction should be refined. The modulated interaction 
obtained here by the proposed one-parameter fitting can be successfully used 
for the Rclativistic Hartree-Bogoliubov calculation. But its applicability to the 
Relativistic Hartree-Fock-Bogoliubov calculation is limited to low densities. 
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Table 1 



Classification of various relativistic approaches to nucleon-nucleon pairing 



Type 


p-h channel 


p-p channel(lowest order) 


References 


B 


DBHF 


bare 


[2] 


B' 


DBHF 


effective range 


[19] 


PI 


RMF 


RMF 


[5-8,23] 


P2a 


RMF 


another phenomenological 


[5,11] 


P2b 


RMF 


bare 


[11,12] 


P2b' 


RMF 


effective range 
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Fig. 1. (a): Curvature of in Eq.(ll) with respect to the cutoff parameter A. (b): 
A-dependence of the pairing gap at the Fermi surface, kp = 0.9 fm~^. These are the 
results of the Hartree-Bogoliubov calculation. 
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Coherence length (Hartree) 
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Fig. 2. (a) Pairing gap at the Fermi surface, and (b) coherence length, as functions 
of the Fermi momentum kp, obtained by adopting the Bonn-B potential, the sudden 
cutoff in Ref . [23] and the two types of effective interaction constructed in this study. 
These are the results of the Hartree-Bogoliubov calculation. 
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Fig. 3. (a) Pair wave function, (b) pairing gap, and (c) matrix element v{kF,k), as 
functions of the momentum k, calculated at a Fermi momentum kp = 0.9 fm~^, 
by adopting the Bonn-B potential, the sudden cutoff in Ref.[23] and the effective 
interaction involving the optimal monopole form factor constructed in this study. 
These are the results of the Hartree-Bogoliubov calculation. 
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Fig. 4. (a) Pair wave function, and (b) pairing gap, as functions of the distance 
r, calculated at a Fermi momentum kp = 0.9 fm~^, by adopting the Bonn-B 
potential, the sudden cutoff in Ref.[23] and the effective interaction involving the 
optimal monopole form factor constructed in this study. These are the results of the 
Hartree-BogoUubov calculation. 
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Fig. 5. Feynman diagram representing the nucleon selfenergy in the Hartree-Fock 
model. 
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Fig. 6. Pairing gap at the Fermi surface as a function of the Fermi momentum kp, 
obtained by the Hartree-Fock-Bogohubov and the Hartree-Bogoliubov calculations 
with a large enough momentum cutoff, 20 fm~^. 
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Fig. 7. The same as Fig. 1(b) but for the Hartree-Fock-Bogoliubov calculations. 
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Fig. 8. Matrix element v{kF,k), as a function of the momentum k, obtained by 
(a) the Hartrcc-Fock, and (b) the Hartrcc-Fock-Bogoliubov calculations with the 
optimal monopole form factors, calculated at three Fermi momenta k-p- 
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Gap at Fermi surface (Hartree-Fock) 
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Fig. 9. The same as Fig.2 but for the Hartree-Fock-BogoHubov calculations. 
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Fig. 10. Energy density including the pairing contributions (both in the 
Hartree-Fock-Bogoliubov and the Hartree-Bogoliubov calculations) and the cutoff 
effects (only in the former) as a function of the Fermi momentum kp. 
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Fig. 11. Energy density at the saturation density corresponding to kp obtained 
by the Hartree-Fock calculations (without pairing), calculated as a function of the 
momentum-cutoff parameter in the form factors, fcp = 1.42 fm~^ for the first two 
cases, while 1.30 fm~^ for the third case. Note that the scale of the abscissa is 
different from Fig. 7. 
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